Tomographic test of Bell's inequality for a time-delocalized single photon 
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Time-domain balanced homodyne detection is performed on two well-separated temporal modes 
sharing a single photon. The reconstructed density matrix of the two-mode system is used to 
prove and quantify its entangled nature, while the Wigner function is employed for an innovative 
tomographic test of Bell's inequality based on the theoretical proposal by Banaszek and Wodkiewicz 
(Phys. Rev. Lett. 82, 2009 (1999)). Provided some auxiliary assumptions are made, a clear violation 
of Banaszek-Bell's inequality is found. 
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I. INTRODUCTION 

The concept of entanglement, as first introduced by 
Schrodinger 1] and discussed by Einstein, Podolsky, and 
Rosen in 1935, has historically been associated with 
systems of two or more quanta. Its counterintuitive con- 
sequences found a mathematical formulation in the work 
of Bell of 1965: pairs of quanta entangled in discrete 
variables (i.e., Bohm-like entangled systems 0) may give 
rise to purely quantum nonlocal correlations. Unfortu- 
nately, no loophole-free test of Bell's inequality has been 
realized so far . Besides their relevance in fundamental 
physics, these phenomena have attracted much attention 
due to their usefulness in quantum information technol- 
ogy. Extravagant but promising protocols such as quan- 
tum teleportation, quantum cryptography, and quantum 
computation have been proposed and experimentally ver- 
ified (see, e.g., and references therein). 

In the early nineties, increasing attention has been 
given to a new perspective of quantum entanglement. 
The concept of entanglement has indeed been extended 
to any system containing a fixed number = 1,2,..) 
of photons, provided at least one pair of spatiotemporal 
modes is involved. Then, even for = 1, one may talk of 
"single-photon n-mode entanglement", provided n > 2. 
The first steps in this direction were taken by Tan, Walls, 
and Collet '7^ (TWC): when a single photon impinges on 
a beam-splitter whose outputs are described by the spa- 
tial modes A and B, the emerging single photon is de- 
scribed by the state \^) — a|l)^|0)B -|-/3|0)a|1)b, where 
|0) is the vacuum, |1) is a single-photon Fock state, and a 
and /3 are complex amplitudes such that |ap -I- = 1. 
In this perspective, a single photon, with its presence or 
absence from a given spatiotemporal mode, may repre- 
sent a tool for entangling two (or more) well-separated 
field modes. In this case, an experimental test of nonlo- 
cality would consist in performing simultaneous measure- 
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ments on the separate field modes "sharing" the single 
photon and check if their mutual correlations violate a 
Bell's inequality. 

Intensive studies and debates have been dedicated to 
the meaningfulness of extending both the concepts of en- 
tanglement and nonlocality to TWC-type s ingle-ph oton 
two-mode systems (see, e.g., i iM O El El 111 I!!). 
Quite recently, Babichev, et al.[ll| and Zavatta, et al.[lj| 
have experimentally characterized a source of two-mode 
spatially and temporally delocalized single photons, re- 
spectively, by employing time-domain homodyne tomog- 
raphy |16l] and reconstructing both density matrix and 
Wigner function of the measured system. In this respect, 
it is worth reminding that the density matrix or Wigner 
function of a given system represents its most complete 
characterization and enables predicting the result of any 
possible measurement one may perform on the system 
]u\ . The results indicate the existence of strong correla- 
tions between the two spatial/temporal modes carrying 
a coherently delocalized single photon. Whether or not 
such correlations are nonlocal cannot be asserted; for in- 
stance, the Bell's inequality test performed in Ref.fTH 
is based upon dichotomization of quadrature data and 
heavily relies upon data-rejection. Of particular inter- 
est, in this respect, is the experiment by Hessmo et al. 
[l3 | , who have implemented the theoretical proposals by 
Tan et al. and by Hardy [3, demonstrating the non- 
classical character of the correlations characterizing two 
spatial modes sharing a single photon. 

Similar to the two-photon case, also "entangled" 
single-photon sources may find immediate application 
in quantum information technology: single-photon ebits 
have indeed been proven to enable linear optics quantum 
teleportation IT^ and play a central role in linear 
optics quantum computation |2(]l l2l]|. Furthermore, the 
time-bin entanglement of the kind of Ref.^T^ (see also 
22, 23J) has been proven suitable for long distance ap- 
plications 13 1^ J where the insensitivity to both depo- 
larization and polarization fluctuations becomes a strong 
requirement. 

In this paper, we wish to make one step further in 
the direction of understanding and proving single pho- 
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ton N-mode entanglement. To this end, we employ the 
single-photon source presented in Ref.jlJ] (see Sectllll|l. 
perform two-mode homodyne tomography (Sect Uvjl and 
employ the reconstructed two-mode density matrix and 
Wigner function to quantify the degree of entanglement 
(Sect. IIV A(l and implement a novel kind of tomographic 
Bell's inequality test (Sects. ITTl and llV B|l . In particu- 
lar, we evaluate the degree of entanglement character- 
izing our single-photon system by employing Vidal and 
Werner's reformulation [2J| of Peres separability criterion 
|26l |. Then, we employ the experimentally reconstructed 
two-mode Wigner function to perform an innovative con- 
tinuous variable Bell's inequality test, based on the theo- 
retical proposal by Banaszek and Wodkiewicz ^3 , which 
we briefly introduce in the next section. 



Our ultra-fast time-domain homodyne detection scheme 
has the potential to achieve this goal psL |2^ . 

It is certainly true that, different from standard Bell's 
inequality tests, our tomographic test imposes an indi- 
rect approach to Bell's inequality: rather than directly 
employing the results of correlation measurements, our 
scheme requires the manipulation of the experimental 
quadrature data in order to reconstruct the Wigner func- 
tion entering Eq. (Q). Note, however, that the reconstruc- 
tion procedure is rather transparent, does not imply ad- 
ditional hypotheses on the state under study, and can 
only introduce noise, without hiding or enhancing any 
information contained in the data. Then, finding a viola- 
tion of Eq. through a tomographic approach indicates 
that also direct measurements would lead to a violation. 



II. TOMOGRAPHIC TEST OF BELL'S 
INEQUALITY: BASIC IDEA 

Let us start by briefly reviewing the theoretical pro- 
posal of Ref . [l0| : the strict connection between the dis- 
placed parity operator (a dichotomic observable) and the 
Wigner function suggests that the Wigner function 
of any two- mode/two-particle system (whether positive 
or negative) may play the role of a nonlocal correlation 
function. Bell's inequality can thus be recast in the form; 

\B\ = —\W{Q,Q)+W{ai,Q)^W{Q,a2)-W{ai,a2)\ < 2, 

. 

where W{ai^ 02) is the value of the two-mode Wigner 
function at the phase-space point (ai = a;i -|- iyi,a2 — 
X2 + iy2) defined by the quadratures Xi and yi of the 
i-ih mode/particle {i — 1,2). This inequality, which 
we shall name Banaszek-Bell's inequality, applies both 
to two-particle and to two-mode entangled systems, in- 
cluding spatially/temporally delocalized single photons. 
Interestingly, compared to all Bell's inequalities theoreti- 
cally proposed for continuous variables and based upon a- 
posteriori dichotomization of the measured field quadra- 
tures |23|, Banaszek-Bell's inequality may achieve higher 
levels of violation. For instance, based on the predic- 
tions of quantum mechanics, the time-bin encoded single- 
photon state [T3 |: 

|*)--i=(|l)„|0)„+i + |0)„|l)„+i), (2) 

where n denotes a well-defined temporal mode (or time 
bin), is expected to maximally violate the inequality of 
Eq. O, giving \B\ w 2.2 for ai = as ~ 0.3. 

The experimental test of Banaszek-Bell's inequality, in 
the form written in Eq. is feasible: our idea is to 
perform two-mode homodyne detection on the two dis- 
tinct time bins (n and n -\- \) carrying the temporally- 
delocalized single photon of Eq. ((2Jl, and to use the ex- 
perimental data to tomographically reconstruct the two- 
mode Wigner function entering the inequality of Eq. 



III. EXPERIMENTAL SETUP 

Let us now consider the experimental setup: a mode- 
locked Ti:sapphire laser emitting 1.5 ps pulses at 786 nm 
at a repetition rate of 82 MHz, is frequency doubled in 
a LBO crystal; the resulting train of phase-locked pulses 
impinges on a nonlinear BBO crystal, cut for degener- 
ate non-collinear type-I SPDC (Sflj]. Signal-idler photon 
pairs centered around 786 nm are generated in symmet- 
ric directions. Before entering an unbalanced fiber-based 
interferometer and being detected by a single-photon 
counter, the idler (trigger) photons undergo a narrow 
spectral and spatial selection aiming to the conditional 
generation of a pure single-photon state in the signal 
mode (see, e.g., [2^ 1^)- Our setup is characterized 
by a single-photon preparation efficiency rjp = .85, which 
depends on both the dark counts of the trigger detector 
and on the purity of the prepared single photon [32| . The 
conditionally-prepared signal beam impinges on a 50-50 
beam splitter (BS) where it is mixed with a strong local 
oscillator (LO) made of an attenuated version of the train 
of laser pulses. Two-mode, high-frequency, time-domain, 
balanced homodyne detection (HD) is then performed 
on two consecutive signal pulses, as depicted in Fig.^). 
The detection efficiency {-qd — 0.74) depends both on 
the efficiency of the two photodiodes (jypD = .88) and 
optical losses, and on the mode-matching with the local 
oscillator {rjMM = .86). The overall experimental effi- 
ciency of our setup is thus expected to be of the order of: 
V = rjpVd = -63. 

From the preparation viewpoint, the key part of our 
setup is the interferometer inserted in the idler chan- 
nel: by setting the interferometer delay T equal to the 
fixed inter-pulse separation Tp characterizing the train 
of phase-locked pump pulses, a click in the idler chan- 
nel does not distinguish the idler photon generated by 
the n*^ pump pulse which travelled the long arm, from 
the idler photon generated by the (n -I- 1)*^ pump pulse 
which travelled the short arm (see Fig. ^)). For an infi- 
nite train of mode-locked pump pulses and in the case of 
equal losses in the two arms of the interferometer, a click 
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a) b) 

FIG. 1: (color online), a) Scheme for the remote preparation 
of a temporally delocalized single photon; i) and ii) indicate 
the two indistinguishable alternatives leading to the coherent 
superposition of Eq. |(2J in the signal arm. b) Scheme for 
the two-mode, time-domain, homodyne tomography of such 
a time-delocalized single-photon state. 

in the trigger detector projects the signal single photon 
onto the coherent superposition of Eq. 

IV. EXPERIMENTAL RESULTS 

In order to check the degree of entanglement and to 
perform the tomographic test of Bell's inequality on 
such conditionally and remotely prepared temporally- 
delocalized single photons, we have performed two- mode 
homodyne detection on the pair of consecutive signal 
time bins sharing the single photon; this is pictorially 
drawn in Fig. ^p) . By making use of what we have re- 
cently named "remote homodyne tomography" to 
remotely vary the relative phase between the two-mode 
state and the LO pulses, we have performed about 10^ 
quadrature measurements on the two consecutive signal 
temporal modes carrying the delocalized single photon. 
This has been made possible by the ultrafast operation 
of the homodyne detector recently developed in our lab- 
oratory ,28, 111- 

The measured field quadratures have been employed to 
reconstruct the density matrix of the measured system by 
means of quantum tomography. In particular, we employ 
the pattern function method (PF) proposed by D'Ariano, 
et al. 33] to directly retrieve the elements of the two- 
mode density matrix in the number-state base from the 
measured quadrature data: pkimn — (fc,/|/5|m, n) is ob- 
tained from the statistical average of the correspond- 
ing pattern functions fkmix,9) over all homodyne data, 
which is: 

Pklmn = {fkni{xi,0l)fln{x2,O2)) (3) 

where xi^2 are the quadratures measured at phases 611^2 
on the two consecutive time bins. The pattern func- 



tions have the form fkm{.x,9) = Fkm{x)e '^"^ with 
Fkm{x) = d {■ipk{x)(f>m{x)) / dx (m > fc), where tpkix) 
and 4>m{x) are, respectively, the regular and i rreg ular 
wave functions of the harmonic oscillator 0, In 
the experiment, 02 — di is remotely varied by means of 
the interferometer while 62 + Oi is left random since our 
system is independent of the global phase. 

When reconstructing the quantum state, we allow each 
single-mode Fock state to contain from zero to two pho- 
tons; each matrix index is thus varied from zero to two 
and a total of 3** = 81 density matrix elements are recon- 
structed for the two-mode system, as shown in Fig. 
No assumptions are made on the system to be recon- 
structed; we just impose a truncation to the reconstruc- 
tion space since no multi-photon contributions are ex- 
pected. 

The density matrix of the measured system has also 
been reconstructed by means of the maximum-likelihood 
method (ML) 35] : the density matrix p that most likely 
represents the homodyne data is retrieved by maximizing 
a functional C{p) involving the POVM associated with 
two-mode homodyne measurements. The only constraint 
we impose is that the density matrix p is a positive- 
definite hermitian matrix with unitary trace. In order 
to limit the number of free parameters in the minimiza- 
tion procedure, we impose p to be independent of the 
global phase 6*2 + 6*1. In principle, the ML method en- 
ables taking into account the experimental imperfections, 
thus reconstructing the density matrix of the measured 
system corrected for the experimental efficiency ry. At 
this stage, we do not impose any such correction. 

Both the PF and the ML method (without correcting 
for 77) give very similar results; hence, in Fig. we just 
plot the elements obtained from the PF method. The 
reconstructed density matrix contains both the expected 
state of Eq. (0) and a vacuum component (poooo) com- 
ing from the non-unitary experimental efficiency. The 
measured system is thus given by the mixture 

= (l-7?)|0),,(0|-f77l*)s.(*|, (4) 

with Y^)s as given in Eq. jSJl; note that almost no multi- 
photon contribution exists in Fig.[2Ji), as expected. From 
the vacuum component of the reconstructed density ma- 
trix we evaluate the overall efficiency to be 77 = 0.61, a re- 
sult in good agreement with our estimated single-photon 
preparation and detection efficiencies. 

A. Degree of entanglement 

The reconstructed density matrix will now be em- 
ployed to quantify the amount of entanglement charac- 
terizing the measured two-mode system. To this end, 
we use the logarithmic negativity Ej^ defined in on 
the line of Peres separability criterion '2^. In general, 
Peres criterion gives a necessary condition for separa- 
bility p^ : if a bipartite system is separable, the par- 
tial transpose of its density matrix (i.e., the transpose 
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FIG. 2: (color online), a) Experimentally reconstructed den- 
sity matrix elements (absolute values) relative to the mixed 
state of Eq. l(4j; b) same elements after "vacuum cleaning" 
by an inverse Bernoulli transformation; c) density matrix el- 
ements reconstructed with the maximum-likelihood method 
incorporating the experimental inefficiency. 



fF^ of subsystem A alone) remains a physical density 
matrix with positive eigenvalues. Then, if we find that 
the partial transpose of our density matrix is not posi- 
tive definite, we know our single-photon system is cer- 
tainly non separable. In order to quantify the degree of 
non-separability (i.e., the degree of entanglement), one 
can evaluate the logarithmic negativity defined as [25l| : 
Ej^{p) = log2 ||yO"^'*||i, where the symbol ||.||i indicates 
the trace norm: ||(t||i = TrV a'^a, which, for Hermitian 
matrices, reduces to Tr\d\ = 1-1- 2\M\, with M sum of 
the negative eigenvalues of a [S^. In other words, the 
trace norm accounts for the negative eigenvalues of fp"'^ , 
and is thus directly related to the amount by which the 
partial transpose fails to be positive definite. In fact, for 
a separable system, ||p'^'^||i = 1 and the corresponding 
logarithmic negativity is zero; for a maximally entangled 
2x2 bipartite system, | \p^^ I |i = 2 and the corresponding 
logarithmic negativity is one. 

In our case, the reconstructed density matrix yields a 
logarithmic negativity £V(/5) — 0.404±0.001, thus show- 
ing the presence of some amount of entanglement even for 
non perfect detection efficiency. This result is perfectly in 
line with the value expected from Eq. with 77 = 0.61. 
Interestingly, the bipartite system described by Eq. is 
one of those "striking" mixtures which remain insepara- 
ble (i.e., fP"^ has a negative eigenvalue) for any value of 
the efficiency rj ^| . A coherently delocalized single pho- 



ton is thus extremely robust against losses: its degree of 
entanglement may decrease but never vanishes. 

B. Test of Banaszek-Bell's inequality 

Based on the above result, our single-photon two-mode 
system is characterized by purely quantum correlations. 
Our next step is to check whether such quantum cor- 
relations are nonlocal. To this end, we reconstruct the 
two- mode Wigner function of the measured system : 

Pklmn 

k,l,rn,n 

where Wij{x,y) is the Wigner function associated with 
the projector \i){j\, and pkimn is the generic element of 
the experimentally reconstructed density matrix. The re- 
constructed Wigner function can now be used to evaluate 
the combination defined in Eq. QJ, thus checking for a 
violation of Banaszek-Bell's inequality. We find that the 
parameter B falls well within the limits imposed by local 
hidden- variable theories (see Fig.|2Jl, a result that may be 
associated with the low overall experimental efficiency. In 
fact, in the case of limited efficiency, the expected Wigner 
function for the state of Eq. takes the form: 

W,ia,,a2) ^4[1 + 2,(K+..P-1)] 

(6) 

hence, taking ai = 0:2 = \fj , the combination given in 
Eq. becomes: 

(7) 

Its behavior is plotted in Fig. O as a function of the 
squared amplitude J for different values of the global 
efficiency rj and it is seen to closely match the experi- 
mental results (squared points) for 77 = 61%. It is evident 
that the current level of experimental efficiency rules out 
the possibility of a loophole-free test of Banaszek-Bell's 
inequality, which would be attainable only for global ex- 
perimental efficiencies larger than 96%. 

To satisfy such strict requirements we need to intro- 
duce some auxiliary assumptions. In this perspective, it 
is worth noting that, different from photon counting, ho- 
modyne tomography allows one to explicitly see the effect 
of the experimental inefficiencies on the measured sys- 
tem: the overall efficiency rj enters into the reconstructed 
density matrix (FiglSt) in the form of vacuum. In par- 
ticular, we find that the vacuum component remains ex- 
actly the same when no interferometer is inserted in the 
idler channel and a single signal photon is perfectly local- 
ized within a given temporal mode. This indicates that 
the non-unitary efficiency t] represents our inability of 
producing and detecting a pure single-photon state, but, 
once such imperfect single photon is coherently delocal- 
ized between separate time-bins, the coherence is not fur- 
ther affected by losses. This is apparent from Fig. |2^), 
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where all the density matrix elements associated with 
the entangled state of Eq. |(2Jl have approximately equal 
weights. In other words, both detection and preparation 
inefficiencies give rise to state-independent losses as if a 
beam splitter with transmissivity 77 were introduced in 
the signal path, mixing vacuum with the original state 
(see Eq. Q and Fig. |21 a)). Correcting for our non- 
unitary efficiencies is thus similar to making the standard 
fair-sampling assumption [^. 

Based on this reasoning, we have chosen to reconstruct 
the density matrix of the measured system by correct- 
ing for the non-unitary efficiency; this is expected to 
remove the vacuum component from the reconstructed 
density matrix, which is, to make the poooo element al- 
most equal to zero and rescale the other elements ac- 
cordingly We will implement such vacuum removal by 
employing two different strategies. Banaszek-Bell's in- 
equality will thus be tested under the fair-sampling as- 
sumption or, more precisely, on the resulting "vacuum- 
cleaned" Wigner function. 

In the first case, we account for losses by making 
an inverse Bernoulli transformation (IBT) of the recon- 
structed density matrix; the result is shown in Fig. [Sis). 
The basic idea is to model the non-perfect experimental 
efficiency by a beam splitter with transmissivity rj fol- 
lowed by an ideal detector: the measured system pmcas 
can be seen as the attenuated version of the system in- 
cident on such fictitious beam splitter. Since the output 
state of the beam splitter is related to the input state by 
a Bernoulli transformation, inversion of this transforma- 
tion gives the state p of the "incident" system, before its 
mixing with vacuum (i.e., before losses affect it) |3^ . 

In the second case, we directly reconstruct the loss- 
free density matrix by adopting the maximum-likelihood 
method (ML) while taking into account the experimental 
efficiency rj — 0.61. The results are shown in Fig.[2t) [37l| . 

As shown in Fig.O the Wigner functions reconstructed 
from both the "vacuum-cleaned" density matrices give 
rise to a good agreement between the experimental data 
and the theoretical prediction of Eq. (0, with unitary ef- 
ficiency. The temporally-delocalized single photon thus 
clearly violates the lower bound of Banaszck-BcU inequal- 
ity under our version of the fair-sampling assumption. 

Despite the observed violation of Banaszek-Bell's in- 
equality, we cannot claim nonlocality for our single- 
photon time-encoded two- mode system. In fact, even 
in the case of unitary efficiency, our measuring scheme 
would still not satisfy the locality hypothesis: quadra- 
ture measurements on the two modes are actually done 
in sequence, so that a physical signal can in principle 
be exchanged between the two co-propagating temporal 
modes while homodyne measurements are performed on 
them. The existence of this possibility makes our Bell's 
inequality test subject to the so-called locality loophole, 
which automatically rules out the possibility of claiming 
non-locality |38l |. However, this loophole simply derives 
from our choice of a temporal delocalization of the single 
photon and can in principle be eliminated in future exper- 
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FIG. 3: Plot of the parameter B,,(J) (Eq. 0) for three 
different values of tlie overall efficiency rj. Continuous curves 
are theoretical predictions, while the points come from the 
experimental data. The horizontal line indicates the lower 
bound imposed by Banaszek-Bell's inequality. 



iments by introducing a fast switch which converts the 
two co-propagating temporal modes into two spatially- 
separated modes; simultaneous measurements could then 
be performed on each of the two time-bins by two distant 
homodyne detectors. 



V. DISCUSSIONS AND CONCLUSIONS 

Our results demonstrate the feasibility of the tomo- 
graphic approach to Bell's inequality and may open the 
way toward new experiments for testing the foundations 
of quantum mechanics. Most important, our experi- 
ment sheds some light on the hotly debated concepts 
of single-photon two-mode entanglement and nonlocality 
0, S M [13, 111 [H El On the one hand, single- 

photon two-mode states of the kind described by Eq. Q 
are known to give rise to interference upon recombina- 
tion of the two spatiotemporal modes through an inter- 
ferometric scheme; in this perspective, the TWC entan- 
glement between two spatial modes sharing a single pho- 
ton can be visualized in terms of wave-particle duality: 
the interference upon recombination supports the wave 
picture of light, while the lack of coincidence counts at 
the two exits of the entrance beam-splitter brings to ev- 
idence its particle nature. On the other hand, we have 
shown that our single-photon two-mode system manifests 
entanglement, both in terms of Peres-like separability 
criteria and, under additional hypotheses, in the more 
restrictive terms of Banaszek's version of Bell's inequal- 
ity. The description of single-photon coherence effects 
in terms of single-photon N-mode entanglement is thus 
certainly possible and meaningful, and could contribute 
to make the counterintuitive wave-particle duality more 
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intuitively accessible. 

In summary, we have evaluated the degree of entangle- 
ment of a single photon coherently delocalized between 
two distinct temporal modes and have implemented an 
innovative test of Bell's inequality based on quantum ho- 
modyne tomography. Wc have shown that, within some 
auxiliary assumptions, the correlations between two well- 
separated temporal modes sharing a single photon violate 
Banaszek-Bell's inequality. In this respect, it is worth 
emphasizing that none of our additional hypotheses is, 
in principle, indispensable. In fact, different from Bell's 
tests relying on post-selection, an improvement in the ex- 
perimental efficiencies would allow us to avoid the fair 
sampling assumption; indeed, a setup with an almost 
perfect single-photon preparation efficiency is currently 
possible. Furthermore a slight modification of the experi- 



mental setup would easily eliminate the locality loophole. 
Thanks to the high efhciencies achievable by homodyne 
detectors, out tomographic approach may represent the 
first step in the direction of a loophole-free test of Bell's 
inequality. In this perspective, also the debated issue 
concerning the possibility of extending both concepts of 
entanglement and non-locality to a coherently delocal- 
ized single photon, finds here an interesting preliminary 
answer. 

This work has been performed in the frame of the 
"Spettroscopia laser e ottica quantistica" project of 
the Physics Department of the University of Florence 
and partially supported by Ente Cassa di Risparmio di 
Firenze and MIUR, under the PRIN initiative and FIRB 
contract RBNE01KZ94. 
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